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Abstract. We study the volume of nodal sets for eigenfunctions of the 
£ — ■ Laplacian on the standard torus in two or more dimensions. We consider 

<^^> a sequence of eigenvalues 4ty 2 E with growing multiplicity J\f — > oo, and 

. compute the expectation and variance of the volume of the nodal set 

with respect to a Gaussian probability measure on the eigenspaces. We 
O |' show that the expected volume of the nodal set is const^/E. Our main 

D ■ result is that the variance of the volume normalized by \fE is bounded 

^0 ' by 0(1/Vf7), so that the normalized volume has vanishing fluctuations 

^-j- , as we increase the dimension of the eigenspace. 

(N 
(N 

1. Introduction 

t-h ; 

The nodal set of a function on a manifold is the set of points where 
it vanishes. Nodal sets for eigenfunctions of the Laplacian on a smooth, 
• compact Riemannian manifold have been studied intensively for some time 

now. For instance, it is known [6 J that except for a subset of lower dimension, 
the nodal sets of eigenfunctions are smooth manifolds of codimension one 
in the ambient manifold. In particular one can define their hypersurface 
volume (in two dimensions this is the length). A conjecture of Yau is that the 
volume of the nodal set is bounded above and below by constant multiples 
of square root of the Laplace eigenvalue. Yau's conjecture was proven for 
real-analytic metrics by Donnelly and Fefferman [7]. The lower bound in the 
case of smooth surfaces is due to Briining [3], see also [3] for planar domains. 

In this paper we study the volume of nodal sets for eigenfunctions of 
the Laplacian on the standard flat torus T d = M d fZ d , d > 2. We write 
the eigenvalue equation as A/ = —47r 2 Ef, where E > is an integer. 
The eigenvalues on the torus always have multiplicities, with the dimension 
J\f = M(E) of an eigenspace corresponding to eigenvalue 4ir 2 E being the 
number of integer vectors A € Z rf so that |A| 2 = E. In dimension d > 5 

this number grows as E — > oo roughly as Ei , but for small values of d, 
particularly for d = 2, the behaviour is more erratic, and depends on the 
prime decomposition of E. 

We will consider random eigenfunctions on the torus, that is random linear 
combinations 

(1.1) f(x) = -i= ^2 6 A cos27r(A,x) - CAsin27r(A,a;) 

^ ^ \& d :\\^=E 
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with b\, c\ ~ N(0, 1) real Gaussians of zero mean and variance 1 which are 
independent save for the relations 6_a = b\, c_\ = — c\. Let £ = Ee be the 
eigenspace associated to the eigenvalue \tx 2 E (i.e. the space of functions of 
form (jl.ip ). We denote by E(«) the expected value of the quantity • in this 
ensemble. For instance, the expected amplitude of / is E(|/(x)| 2 ) = 1. 

Denote by Z(f) the volume of the nodal set of an eigenfunction (11.11) . 
Our first result, Proposition 14. 1\ is that the expected value of Z is 

E(Z) = const ■ VE 

for a certain constant depending only on the dimension d. This is of course 
consistent with the bounds of Donnelly and Fefferman [7] . 

Our main result, Theorem 16. 1\ is that the variance of the normalized 
volume Z/ \[E is bounded by 

Z 1 
Var(^=) <C —=, as N — ► oo . 
WE 

(We believe that the correct upper bound for the variance is 0(1/ N)). 
Thus the fluctuations of Z(f)/y/E around its mean value die out as the 
multiplicity J\f tends to oo. Note however that Z(f)/yE is not asymp- 
totically constant; for instance, if E = dm 2 then for the eigenfunction 
f( x ) = rij=i sm ^ 7rmx i we have Z(f)/\/ r E = 2y/d while if E = m 2 then 
for the eigenfunction f(x) = sin2-7rma;i we have Z(f)/y/~E = 2. 

Theorem l6.1l can be viewed as lending support to the expectation^! that for 
eigenfunctions on negatively curved manifolds, which are believed to behave 
similarly to random waves [2], the volumes of nodal sets, normalized by the 
square-root of the eigenvalue, do tend to a limiting value. See [T2J for some 
work on the complexified nodal set of eigenfunctions in this context. 

Previous work in this vein is due to Berard [I], who computed the ex- 
pected surface measure of the nodal set for eigenfunctions of the Laplacian 
on spheres. Neuheisel |10j also worked on the sphere and gave an upper 
bound for the variance. Berry [3] computed the expected length of nodal 
lines for isotropic, monochromatic random waves in the plane, which are 
eigenfunctions of the Laplacian with eigenvalue E. He found that the ex- 
pected length (per unit area) is again of size about \[E and argued that the 
variance should be of order log E. 

More recently, F. Oravecz and ourselves have investigated a different char- 
acteristic of the nodal set of eigenfunctions on the torus, namely the Leray 
nodal measure [11] , and have succeeded in obtaining the precise asymptotics 
of the variance of the Leray measure as M — > oo. 

1.1. Plan of the paper. We employ a version of the Kac-Rice formula for 
the volume of the nodal set, which using the Dirac delta function can be 
written as 

Z(f) = I 8(f(x))\Vf(x)\dx , 

see $3]for the rigorous version. To compute the expected value of Z is then a 
simple matter once we find that f(x) and df /dxj are independent Gaussians. 
This is done in § UJ In § [5] we derive a formula for the second moment of Z, 



We thank Steve Zelditch for a discussion of this. 
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which requires knowing the covariance structure of the 2d + 2-dimensional 
Gaussian vector v(x,y) = (f(x), f(y), V/(x), V/(y)). That v(x,y) is indeed 
a non-degenerate 2d + 2 dimensional Gaussian is verified in the Appendix. 
As a result, we find that E(Z 2 ) = fj d K(z)dz, with 

1 f exp(— \v Q(z)~ 1 v T ) dv 



K(z) = — = / IfiI]^ 



v/1 " ^) 2 ■/*" y/detQ(z) (27T) d + 1 ' 

where u(z) = K(f(x)f(x + z)) is the two-point function of the ensemble, and 
where J)(z) is a certain positive definite 2d x 2cf matrix which enters into 
the covariance structure of the Gaussian vector v(x,y). In section 53 which 
is the heart of the paper, we bound the variance of Z. 

2. The model: random eigenfunctions on the torus 

2.1. Random eigenfunctions. We consider non-constant eigenfunctions 
of the Laplacian on the standard flat torus T d = M^/Z^. The solutions of 
the eigenvalue equation 

Aip + 4ir 2 E?p = , E ^ , 

form a finite dimensional vector space £ = having as a basis the expo- 
nentials e 2m ( X ' x } , for A in the frequency set 

A = A E = {\eZ d ,\X\ 2 = E} . 

We define an ensemble of Gaussian random functions / G £ by 

f(x) = —j= b\ cos 27r(A, x) — c\ sin 27r(A, x) 
^ xTa 

with &A) c a ~ A(0, 1) rea l Gaussians of zero mean and variance 1 which 
are independent save for the relations b_\ = b\, c_\ = —c\. Thus we can 
rewrite 

(2.1) f(x) = J— ^ &a cos 27r(A, x) - c A sin27r(A,2;) 

AeA/± 

where now only independent random variables appear. With our normal- 
ization, we have E(|/(x)| 2 ) = 1 for all x G T d . 

Definition 2.1. An eigenfunction f E £ is sing ular if3x G T d with f(x) = 
and V/(x) = 0. An eigenfunction f G £ is nonsingular if V/ / on t/ie 
nodal set. 

Lemma 2.2 ([llj, Lemma 2.3). The set of singular eigenfunctions has codi- 
mension at least 1 in £ , and so has measure zero in £. 

2.2. Properties of the frequency set. The dimension M = dim£ is the 
number of the frequencies in A, which is the number of ways of expressing 
E as a sum of d integer squares. For d > 5 this grows roughly as E d l 2 ~ 1 
as E —* oo. For d < 4 the dimension of the eigenspace need not grow with 
E. For instance, for d = 2, A/ - is given in terms of the prime decomposition 
of E as follows: If E = 2 a J| • J| fc q^ 7fe where = 1 mod 4 and % = 3 
mod 4 are odd primes, a,(3j,jk > are integers, then N = 4]X(/3j + 1), and 
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otherwise E is not a sum of two squares and N = 0. On average (over 
integers which are sums of two squares) the dimension is const ■ y/logE. 

The frequency set A is invariant under the group Wd of signed permuta- 
tions, consisting of coordinate permutations and sign-change of any coordi- 
nate, e.g. (Ai,A2) •— > (— Ai,A2) (for d = 2). In particular A is symmetric 
under A i — ^ — A and since ^ A, we find M is even. We write A/db to denote 
representatives of the equivalence class of A under A i— ► —A. 

We will need some simple properties of A: 

Lemma 2.3. For any subset O C A which is invariant under the group Wd, 
we have 

( 2 - 2 ) |0| X] A J' Afc = 7 ' 5 i' k ■ 

1 1 xeo 

Moreover for any C G M. d , 

1 1 xeo ' 

Proof. For i ^ j use the symmetry of O under the sign change of the i-th 
coordinate to change variables and deduce that the LHS of (|2.2p vanishes. 
For i = j note that the sum J2xeO A ? * s independent of i since O is symmetric 
under permutations; hence we may average the RHS over i to find that 

xeo i=i ago xeo 

proving (|2.2p . To prove (|2.3[) we expand (C, A) 2 = Yl'j k=i c j c k^j^k and use 
(12T2D . ' □ 

Note that (|2.3I) implies that the frequency set A spans M d . 

2.3. The two point function. The two-point function of the ensemble is 

(2.4) u{z) := E(f(x + z)f(x)) = A ^ cos 2tt(A ! z) . 

AeA/± 

The two-point function clearly satisfies \u(z)\ < 1. We will need to know 
some of its basic properties, proved in [TT], which we summarize as: 

Proposition 2.4. The two point function satisfies 

(1) There are only finitely many points x 6 T d where u(x) = ±1. 

(2) The mean square of u is Jj d u 2 = jj- . 

(3) The mean fourth power of u is bounded 6x0 Jjd u A <C 1/Af ■ 

(4) The kernel — u 2 is integrable on T d . 

Part [1] follows from lemma 2.2], part [3] is Proposition 7.1], and 
part [J] is [TH Lemma 5.3]. 



2 Except possibly in dimensions d = 3,4 we have a better bound in [TT] of o(l/Af), 
though we have no use for this finer information in this paper. 
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3. A FORMULA FOR THE VOLUME OF THE NODAL SET 

Let x be the indicator function of the interval [—1,1]. We define for e > 



Lemma 3.1. Suppose that f 6 £ is non-singular. Then 

vol(/- 1 (0)) = limZ e (/). 

e— >0 

Proof. By the co-area formula [8], for / smooth and <fi integrable, we have 



(j)(x)\Vf(x)\dx = / <j)(x)dx ds . 

Taking (j)(x) := tj-x(^~^)> which is constant on the level sets f~ 1 (s) gives 



Uf) = \- e fjol(f-\s))ds. 



Now if / is non-singular then s t—* vol(/ is continuous at s = and so 

by the fundamental theorem of calculus, 

lim 1 ^ vol(f-\s))ds = vol(rHO)) . 

Thus lim^o Z e (f) = vol(/ _1 (0)) as claimed. □ 
Lemma 3.2. For all f € £ we have 

Z e (f) < 6d^/~E . 

We begin with the one variable case which we state as a separate lemma 
(cf Lemma 2]): 

Lemma 3.3. Let g(t) be a trigonometric polynomial of degree at most M . 
Then for all e > we have 

1 



2e 



\g'(t)\dt < 6M . 

{t:|3(*)l<4 



Proof. We partition the set {t : \g(t)\ < e} C [0, 1] into a union of maximal 
closed intervals [afe,6fe] (with a k < bk), disjoint except perhaps for common 
edges, such that on each such interval g 1 has constant sign, that is either 
g' > or g' < 0. If g' > on [a k , bk] then either g(ak) = — e or g'(ak) = 
and a& is a local minimum for g, and < +e. If </ < on [0^,6^] 

then either g(a k ) = +e or g'(a k ) = and a k is a local maximum for g, and 

ffOfc) > - £ - 

If 5' > on [0^, 6^] then 

|5 / (t)|dt= / g'(t)dt = g(b k )-g(a k )<2e, 

while if </ < on [a&, &&] then 

* |^(t)|rft = / " -g{t)dt = g(a k ) - g(b k ) < 2e . 
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Thus the total integral is bounded by the number v of intervals [a^, 

-f 

26 J{t:\g(t)\<e} 



\g'(t)\dt < v . 



Now the number of intervals is bounded by the number of a's for which 
g(a) = ±e plus the number of a's for which g'(a) = 0. Since both g and g' 
are trigonometric polynomials of degree < M, the number of such intervals 
is therefore 3 • 2M = 6M. This gives the required bound. □ 



We now prove Lemma 13.21 by reduction to the one-dimensional case. 
Proof. Since |V/| < Y^j=i l&jfrl we have 



.-i e oxj 



and we will bound each term. Taking j = 1, we have 



Of 



dx 




{teTl : [/(i,jf)[<e} 



df(t,y) 



dt 



dt ) dy . 



In the inner integral we have for each y G T^" 1 a one variable polynomial 
g(t) = f(t,y) of degree at most \[E and hence by Lemma [3.31 the inner 
integral is at most 6\^E. Summing over j introduces another factor of d. □ 

As a consequence of the fact that for nonsingular functions we can com- 
pute the volume Z(f) of the nodal set of / € £ via Lemma l3.1l and the fact 
that almost all / £ £ are nonsingular (Lemma l2.2p . we find: 

Corollary 3.4. The first and second moments of the volume Z(f) of the 
nodal set of f are given by 



E(Z) =E(limZ e ), 

€—*0 



E(Z 2 ) =E( lim Z tl Z e2 ) 

£1,62— »0 



4. The expected volume of the nodal set 
In this section we show 
Proposition 4.1. For d > 1, 

E(Z) = l d VE 

where 

-Lr! — • 



d r(|) 



Proof. Since Z e is uniformly bounded by Lemma 13.21 we can use the Domi- 
nated Convergence Theorem to write 

E(Z) = E(limZ £ ) = limE(Z £ ) . 



VOLUME OF NODAL SETS 



7 



By Fubini's theorem, 

E(Z e ) =^(j e J Td x(^)|V/(*)|dz) 

= / E(± X (—)\Vf(x)\)dx=: f K e {x)dx. 

Now for each x € T d , the function / ^x(^^)|V/(a;)| is bounded and 
hence we may evaluate the integral by using the joint probability density 
of the variable (f(x), V/(x)), whose components are Gaussian of zero mean 
with covariance 

E(/(x) 2 ) = l, E (/(*)^(*>) =° 

and 



(41) E( ^-^r } = ^ AjAfe = 4 " 7 • 

AeA/± 

by Thus 

1 f a 2/ 2 da /" r d|6| 2 d& 

Kefc) = 7T / x(-)e /2 ^= / Hexp(- 



sfd- (2vr)( rf + 1 )/ 2 jR* eXP[ 2 m /u, "2e 
Integrating over T d and taking the limit e — ► gives 

E(Z) =I d VE 

where 

If 1 

= — = ; — / Ivlexpf — \v\ 2 )dv . 

\/d(27r)( d -l)/ 2 Ad ' V 2 1 1 ' 

In the one-dimensional case, Z x = / R |v|e- 1,2 / 2 d?; = 2. For d > 2, 

1 /*°° 

|u| exp(--|£| 2 )d?; = vol{S d - 1 ) / re^/V^dr . 

2 In 



Using 



vol(S d -l) = ^l, r d e-^dr = 2^r(^) 

1 ' T(IY Jo ^ 2 ; 



2 

gives 



[ \v\eM-hv\ 2 )dv = V2(2n) d ^^-l 



(which is consistent with the computation for d = 1). Thus 



1h 



d r(f) 



as claimed. □ 



s 
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5. An integral formula for the second moment 

5.1. The covariance matrix. The covariance matrix y) of the Gauss- 
ian vector {f(x),f(y),Vf(x),Vf(y)) is given by 



(5.1) E = , nT 



where 



A B 
B T C 



_ ( E(/(x) 2 ) E(f(x)f(y))\ _ (E(f(x)Vf(x)) E(/(x)V/(y))\ 
~ \E(f(y)f(x)) E(/(y) 2 ) )>"- \E(f(y)Vf(x) E(f(y)Vf(y))J 

and 

(E(Vf(x) T Vf(x)) E(Vf(xfVf(y))\ 
\E(Vf(yfVf(x)) E(Vf( y fVf(y))) ' 

For generic (x,y), the covariance matrix is nonsingular (see Appen- 

dix|H). 



C 



Lemma 5.1. The covariance matrix T,(x,y) depends only on the difference 
z = x — y and is given in terms of the two-point function u by 

A(z) B( z y 

B(zf C{z) t 



^(x,y) = \ r >,.:\T 



where 



A(z)-[ 1 u{z) ] B(z)-( 6 ~ V ^ z) 



(here 0,S/u are row vectors), and 

C(z) = ( 



^1 -H(z) 



-H(z) 



o2 

where H = ( Q^.g Xk ) is the Hessian of u. 

Proof. By definition of the two point function, we have A = ( ) . To 



compute B, use 



8/ MW ,, d 

and hence 



n^r(x)f(y)) = —E(f(x)f(y)) = d 3 u{x - y) 



dxj dxj 



Hf( x )djf(y)) = dju(y - x) = -dju(x - y) 

In particular 

E(f(x)Vf(x)) = 

Therefore 

B(z)=( 3 

(where denotes the d-dimensional zero row vector) . 
To compute C, use (|4.1|) to find 

E(Vf(xfVf(x)) = ^I d . 
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More generally 



^2 C^^TA 

and so 

E(V/(x) T V/(y)) = "(^^(^ - = -H(x - y) 

Thus 



as claimed. 

The inverse of £ (when it exists) is given by 



s- 1 



* jr 1 



with f2 being the 2d x 2cf matrix 

n = c - B T A~ l B . 

We will call the reduced covariance matrix. We have 

(5.2) detS = detyldetO = (1 -u 2 )detO . 

By Lemma l5.1( we have 

Q _(^ 2 (E/d)I -H \ 1 / D T D uD T D\ 
[b - 6) U ~\ -H 4ir 2 (E/d)l) \-u 2 \uD T D D T D) 

where D(z) = Vii(z) and H = ) is the Hessian of u. 

5.2. A formula for the second moment. 
Proposition 5.2. The second moment of Z(f) is given by 

(5.4) E(Z 2 ) = f K(x)dx 

Jf d 

where 

/, ,x s 1 [ exp(-i?;0- 1 w T ) dv 

(5.5) K (x) = —== / Fi||p 2 " 



v 7 !^?^ 11 11111 VdetQ (27r) d +! • 

Denote 

K eue2 (x,y) : =^^||V/(x)||||V/(y)||x(4f)x(^)^(/) 

We have the following 
Lemma 5.3. For (x, y) £ T d x T d with u(x - yf / 1 

(5.6) K eu62 (x, y) < d 



y^l - u 2 (x - y)' 
where the implied constant depends only on the dimension d. 
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Proof. Write f(x) = (f, U(x)), where U(x) is the unit vector 

a/2 / „ . \ _ , v v.-i 



U(x) = cos 2tt(\, x), sin 2ir(\, x) G S~ 

VAAV / AeA/± 

and where we identify the function / with a vector in R-^ via (pET]) . Note 
that (U(x),U(y)) = u(x — y) is the cosine of the angle between C7(x) and 

We have 

Vf(x) = DU-f 
where the derivative DU is a d x M matrix. Equivalently, 



1< t < d 



By the triangle and Cauchy-Schwartz inequalities, 



||V/(z)||<£ 



i=i 



by a computation of Therefore 



dx, 



(5.7) W* lS )<^ L l<ei ll/H V " /l|2/2 ^ • 

\f{y)\<e 2 

Consider the plane ir C M - ^ spanned by C/(x) and U(y). The domain of 
the integration is all the vectors / E so that the projection of / on tt 
falls into the parallelogram P of lengths 2ei and 2e2- The cosine of the angle 
between the sides of P is (U(x), U(y)) = u(x — y). Therefore the area of P 
is 

area(-P) = 4eie2— ^ . 

1/1 - u(a; - y) 2 

Write the multiple integral in (|5.7I) as the iterated integral 

(5.8) / (7 \\ffe-W fr/2 df)d P , 

where the variable p runs over all the points of the parallelepiped P. The 
inner integral in (|5.8|) is 0(1) with the constant depending on d only. Indeed, 
note that for every /i G tt 1 , 

||p + /l|| 2 e~ l|p+/l||2/2 = (|bf + ll/iU^e^lWP+WI 2 )/ 2 

«(l + ll/i|| 2 )- e HI/l||2/2 , 

since ||p|| 2 e~" p " 2 / 2 is bounded. Our claim follows from convergence of the 
integral J* K Af-2(l + ||HI 2 ) e ~" u '" ^ 2 dw. Therefore 

L)i<« ii^n 2eHI/l|2/2 ^ « area ^ « ^ 7i ! , 2 • 

Substituting the last estimate into (|5.7p proves (|5.6j) . □ 
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Proof of proposition \5.SX By Corollary 13.41 we have 



E(Z 2 



[ lim / 



[|V/(x)|| X (^)dx^- 



\\Vf(y)\\ X (^-)dy 

jd €2 



duU) 



where \i is the Gaussian measure d^,{f) = e We wish to 

change the order of the limit and the integration. To do so, we notice that 
by Lemma [3. 2 1 the integrand is bounded by 0{E). Therefore, the change of 
order follows from the dominated convergence theorem. Thus the integral 
equals 

' ' ' l|V/(.*0||||V/(/;)||\( 

'(T d ) 2 



lim 



ei,€2— »0 4ei£2 j£ , 

Using Fubini's theorem, this equals to 
(5.9) 



tt*)\ (ttv) 



ei 



X 



(2 



dxdydfi(f) . 



1 



lim 



l|V/(x)||||V/(y)||x 



f(*)\ ff(y) 



ei 



X 



(■2 



dfi(f)dxdy . 



'(T d ) 2 J£ 

Now we wish to exchange the order of taking limit and the integration 
over (T d ) 2 . To justify it, we use the dominated convergence theorem with 
Lemma 15.31 The upper bound for u(x) ^ ±u{y) is sufficient, since this 
happens for almost all (x,y) 6 {T d ) 2 , and changing the values of a function 
on a set of measure does not have any impact on the integrability and the 
value of the integral of a function. The convergence of the RHS of (15. 6\\ was 
shown in [IT] . Therefore, we may exchange the order of the limit and the 
integral in (|5.9p to obtain 



(5.10) 
where 



E(Z 2 



K(x,y) 



K(x,y)dxdy, 



T d xT d 



lim K ei , e2 (x,y) . 

61,62— >0 



We will replace the vector (f(x), f(y), Vf(x), V/(y)) with a 2<i+2 dimen- 
sional Gaussian vector with covariance matrix S(x, y) = S(x — y) defined in 
(|5.1|) . The proof that for almost all x, y this is indeed a 2d + 2 dimensional 
process, is relegated to Proposition IA. II in the Appendix. This gives 

(5.11) 

1 f (W\ ^Y w 2^ c - faw )E-ifa w )T /2 dvdw 

v / deTS(27r) d + 1 
dvdw 



4eie 2 
i r e i r e 2 



Vl\\\\V2\\X 



ei J\ £2 

{v,w)j:- 1 (v,w) T /2_ 



4eie 2 j- ei j-e 2 
We therefore have 



pi »2 e 



V / detS(27r) d + 1 



K(x,y) 



lim 



1 



\vi\\ \\V2\\e 



-(v,w)T,- 1 (v,w) T /2_ 



dvdw 



V / detS(2vr) 



d+l 



12 
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Since the last integrand is continuous, we may use the fundamental theorem 
of the calculus to replace the averaging over w\, W2 by the value at w\ = 
W2 = 0, to obtain 

dvdw 



K[x,y) = / |H||| W2 ||e-M)£->,oT/2. 

' - " v / detS(2^) rf + 1 



We have 



s- 1 



* * 

* n- 1 



with Q = C — B T A 1 B being the reduced covariance matrix, which we 
computed in (15. 3p . Together with (|5.2[) we find 

1 f , ex P( _ \vQ.~ 1 v T ) dv 

KM = 7T^L NH Vditn (2^ • 

Finally, we get Proposition 15.21 by noticing that K(x,y) = K{x — y), 
and therefore the (double) integral in (|5.10p may be expressed as a single 
integral. □ 

In the course of the proof we saw that K(x,y) = lim £l)£2 _ > o Kti,e 2 ( x ) ■ 
Therefore, taking the limit ei, e% ^0 and using Lemma 15.31 we obtain 

Corollary 5.4. If u(x) 2 / 1 then 

E 



K(x) < 



y/l — u(x) 2 
6. A BOUND FOR THE VARIANCE 

In this section we prove: 
Theorem 6.1. For d > 2, 

Var(Z) = O ' L 



6.1. Plan of the proof. We use the integral formula (|5.4|) for the second 
moment of Z(f), that is K(Z 2 ) = J T<J K(z)dz, with 



K(z) = — / [|ui||||v2 



exp(— 7}vQ l v T ) dv 
Vl^ 1 "^ lll|iyz|1 ^/ditO (27r) d + 1 ' 

As in [llj . we will define a notion of "singular points" in where the 
factor l/vi — v? is large, and treat separately the singular and nonsingular 
points. The singular set is shown to give a contribution of 0{E/M). On the 
nonsingular set, the factor — u 2 may, up to an admissible error, safely 

be replaced by 1. To treat the Gaussian integral, we write 

47T 2 £' 

n(*) = — (i-s(z)) 

and recover the square of the expected value K(Z) 2 from the contribution 
of the identity matrix I; the rest is then the key quantity for bounding 
the variance. Setting o~(z) to be the spectral norm of S(z), we show that 
that variance is bounded by E (J* Td a(z)dz + 0(l/jV)). Now a(z) is at most 
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\/tr(S(z) 2 ), whose integral we need to bound. We do this by using Cauchy- 
Schwartz, which allows us to bound it by (J Jd tv{S{z) 2 )dz) 1 / 2 <C l/y/J7. 
Hence the variance is Var(Z) <C E/yfJJ. It should be possible to improve 
this to 0(E/N). 

6.2. The singular set. We give the definition of [11] for singular points: 

Definition 6.2. A point x £ T d is a positive singular point if there is a set 
of frequencies A x C A with density Jj^ji > \ _ 1^ f or which cos 2ir(\, x) > 3/4 
for all AgA x . Similarly we define a negative singular point to be a point x 
where there is a set A x C A of density > 1 — ^ for which cos 2ir( A, x) < —3/4 
for all A E A x . 

Let M m \[E be a large integer. We decompose the torus T d = M d /Z d 
as a disjoint union (with boundary overlaps) of M d closed cubes of side 

length 1/M centered at k/M, k E Z d . 

Definition 6.3. A cube I? is a positive (resp. negative) singular cube if it 
contains a positive (resp. negative) singular point. 

Definition 6.4. The singular set B is the union of all singular cubes. 

In [TT], we showed that the measure of the singular set is bounded by 

(6.1) meas(B) < / u(x) 4 dx < —. 

(and except in dimensions d = 3, 4 this is o(l/jV)). 

In order to bound the contribution of the singular set to the integral in 
(|5.4p . we use Corollary 15.41 It was shown in [11] (see (6.3)) that 



(6.2) f dX « I u{xfdx « 1 . 

Jb a/1 - u(x) 2 JT d A/ 



Therefore we obtain: 



Corollary 6.5. The contribution of the singular set is bounded by 

E 



K(x)dx <C El u(x) 4 dx <C 

JT d 



6.3. The nonsingular set. We now want to estimate the contribution of 
the nonsingular set to the integral formula of Proposition 15. 21 for the second 
moment of Z. Recall that it reads K(Z 2 ) = J Td K(z)dz with the kernel K(z) 
given by (15. 5|) . that is 

w x 1 f exp(-i^(z)- 1 t; T ) dv 



v/1 - «(^) 2 11 " x 11 11 " z 11 ^det^(z) (2TT) d + 1 ' 

A consequence of the definition of singular points is that on the nonsin- 
gular set, u is bounded away from ±1. In [TTl lemma 6.5] we showed that if 
x E T d is nonsingular then 

Kx)|<i-^. 
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As a consequence, on the nonsingular set, we may expand 

1 



1 + 0(u 2 



where the implied constant depends only on d. 

We now wish to handle the "reduced covariance matrix" f2 of (j5.3H on 
the nonsingular set. We write O = (Ai: 2 E/d) • ill and tti = I — S, where 

ffi*\ a- d 1 ( ° TD (1-u 2 )H + uD t D 

V- 6 ) b ~ A7T 2 E ! _ u 2 [ {1 _ U 2 )H + uD T D D T D 

Note that since outside a set of measure zero, Q.\ 3> is positive definite, we 
have S <C I in the sense that all eigenvalues of S are in (— oo, 1). Let a be 
the spectral norm of 5, so that denoting the eigenvalues of S by cti, . . . , a2d, 



l<7<2d 

We give a bound on the mean and the mean-square of a on the complement 
B c of the singular set. 

Lemma 6.6. 

(6.4) j B / dx «h 



and 

(6.5) / adx < 



77 



Proof. The bound fl6.5f> follows from (16. 4p by applying Cauchy-Schwartz, so 
it suffices to prove (IfHlh We have a 2 < £ «| = tr(S 2 ), and so it suffices to 
show 

(6.6) / tr(5' 2 )dx <C . 

On the nonsingular set, the expression ^ a is bounded, and hence for 
purposes of upper bounds may be ignored. The entries of S 2 on the nonsin- 
gular set are thus bounded by sums of the following expressions : 

1 d 2 u d 2 u 1 d 2 u du du 1 du du du du 



E 2 dxidxj dxkdxi ' E 2 dxidxj dxk dxi ' E 2 dxi dxj dx^ dx# 

and it suffices to show that the integral of each over all of T d is 0(1/Af). 
By applying Cauchy-Schwartz, it suffices to show 

. 8 2 u , 2j E 2 f du . 4j E 2 

bir ) dx < —r and / — fdx < . 

We have 

9 2 « -Sir 2 



' ** icA/4- 



AeA/± 
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and hence 

}2„, o_2 



A„ueA/± 
AeA 



since A 2 < |A| 2 = E 

f ( dr 

h d \ dx k 

8u — 47T 



To bound Jjd(-§^Y dx , we write 



Afc sin27r(A, x) 



dxk M 

AeA/± 

and as above, we have 

E 2 

A\...,A 4 eA/± 
A 1 ±A 2 ±A 3 ±A 4 =0 



/ (l^) 4 ^ « ^4 E 

7 T /9x fc ; AA 4 ^ 



since A 1 , A 2 , A 3 determine A 4 once we decree that A 1 ± A 2 ± A 3 ± A 4 = 0, and 

\X h \ <^Ve. □ 

6.4. Concluding the proof of Theorem 16.11 Since fii = O • d/(&ir 2 E) 
is symmetric and positive definite (away from a set of measure zero), it has 
a positive definite square root Pi = Pf 3> 0, Oi = P 2 . By proposition 15.21 

zrt \ 1 f ii- .expC-i^" 1 ^) dv 

K \ x ) = n o / \ v l\\ v 2\ 



(zPl)i (zPi) 2 |e |zl /2 



on using the change of variables v = 2n J^ - z - P\. 
We claim that 

Pi = 7(1 + 0(a)) . 

Indeed, if S = UDU T , U orthogonal and D = diag(ai, . . . , a2d) then Pi = 
U(I — D) l l 2 U T and using the inequality \y/\ — a— 1| < \a\ for — oo < a < 1, 
gives 



(I-D) 1 ' 2 = I + 0( 



/l"i| 



\ |«2d| 



7(1 + 0(a)). 



Thus we may write zP\ = z (1 + O(cr)). 
On the nonsingular set, we may expand 

' 1 + 0(u 2 



u 



2 
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and so we find that on the nonsingular set 

= E(Z) 2 (l + 0(u 2 ) + 0(a) + 0(a 2 )) . 
Integrating over the nonsingular set, and using 

/ 1 = 1+ 0(meas(5)) 
JB C 



we find 

f B c K(x)dx 
E(Z 2 ) 



1 + (/ ^ + ° + °" 2 ) cfe ) + °( meas ( B )) 



Now J Td u(x) 2 dx = jj , and by lemma 

7 2 dx <C -77 , / ddx <C — — 
M' J B c v(AT 

meas(-B) <C / ii 4 dx <C v? ; 

Jtd N 



Furthermore, by ()6.1 



so that we find 



K{x)dx = E(Z 2 ) ( 1 + 0(-J=) 



By Corollary 16.51 the singular set contributes at most 

j B K{x)dx«^. 



Therefore we find 

E(Z 2 )=E(Z) 2 + 0(^=), 
VN 

that is 

Var(Z) < L 



IN 

Thus we have concluded the proof of Theorem 16.11 □ 

Appendix A. The non-degeneracy of the covariance matrix 

In this appendix we show that the covariance matrix defined by (|4.1|) is 
nonsingular for almost all (x,y) G (T rf ) 2 , thereby justifying the change of 
variables (|5.1ip . 

Proposition A.l. Assume that N ~^>d 1 an d d > 2. Then for almost all 
(x, y) G T d x T d the linear map £ — ► M. 2d+2 defined by 

f^(f(x),f(y),Vf(x),Vf(y)), 

is surjective. 
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We want to show that for almost all pairs (x, y) G T d x T d , the only vector 
(a, 0, C, D) G R 2d+2 satisfying 

af(x) + Pf(v) + 7^(C,Vf(x)) + ±-(D,Vf(y)) = 0, V/ G £ 
is the zero vector. Taking f(x) = e 2m ^' x \ A G A gives 

ae 2*i<A,x) + p£«i{Kv) + ie 2ni(X,x) ^ ^ + ^(A.J/) ^ A) = 0, VA G A 

or setting z = y — x, 

a + i(C, A) = - e 2?n < A ' 2 > (/? + A)) , VA G A . 
Thus we are reduced to proving the following: 

Lemma A. 2. Assume that J\f 1 and d > 2. Then for almost all z, the 
only solution for the equation 

(A.l) a + i(C,X) = -e 2ni{x ' z) ((3 + i{D,X)), VA G A, 

is a = (3 = 0, C = D = 0. 

Proof. We divide the work into two steps: In the first step, we show that 
for all z G T d , the solutions of (jA.ip satisfy j3 = ±a and Z? = dbC. In the 
second step, we take (3 = dba and Z) = ±C and show that for almost all 
z G T d , the only solutions of (|A. 1[) are a = and C = 0. 

Step 1: We first show that for all z G T d , all solutions of (jA.ip satisfy 
[3 = ±a and D = ±C. 

Taking squared norms of both sides of (|A. 1[) , we get 

a 2 + (C, A) 2 = /3 2 + (£>, A) 2 , VA G A 

or 

(C, A) 2 - (D, A) 2 = /3 2 - a 2 , VA G A . 

Setting A = C — D = (ai, «2, • • • ) and B = C + D = (b\, 62, • • • ), we have 

(A.2) (yl, A) • (ZJ, A) = (3 2 - a 2 , VA G A 

and it suffices to see that A = or B = 0. 

If 1) then there is some A G A with two nonzero coordinates, say 

A1A2 7^ (by applying a permutation of the coordinates to A we may replace 
1 and 2 by any pair of distinct indices). For each e G {±l} rf , replace A in 
(HDby 

X e := (eiAi, 62A2, . . . , e^Arf) , 

multiply the result by 

Xl,2(e) = eie 2 

and sum the resulting equalities over all e G {±l} d , using 

J2 Xi,2(e) = 

eS{±l} d 

to get 

J2 Xi,2(e)<A,A e )-<S,A £ )=0. 
ee{±i} d 
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Expanding 

d 

(A,X e ) ■ (B,X e ) = ^ ejekQjbkXjXk 
j,k=i 



and using 



we get 



Xi,2(e)e i e fc = I 
ee{±i} d I 



2 d , (j,k) = (1,2) or (2,1) 
otherwise 



2 d AiA 2 (ai6 2 + a 2 6i) = 
and since we assume AiA 2 7^ 0, we find ai& 2 + ai6 2 = 0. Repeating the 
argument with any pair of distinct indices finally shows that 

(A. 3) ciibj + djbi = 0, \/i ^ j . 

If A 7^ 0, say ai 7^ 0, then we find that 

(A.4) bj = -^ aj , Vj + 1 . 

ai 

Thus if b\ = then all bj = 0, that is B = and we are done. Therefore we 
may assume that b\ 7^ (and we have also assumed a± 7^ 0). We will show 
this cannot happen. 

If d > 2, we substitute (IA.4j) in (|A.3j) with any i 7^ 1, j 7^ 1 to get 

&i 

2 — a^a,- = 
ai 

that is since 61 7^ 0, that 

didj = 0, Vi^j, 

Thus there is at most one index k 7^ 1 with 7^ 0, say fc = 2, so we find 
that Oj = for j / 1, 2, and by (TOl) we therefore have bj = for j 7^ 1,2. 
Thus 

A=(ai,a 2 ,0), B = — (01,-02,0) 
ai 

(if d = 2 this still holds, we just ignore the extra coordinates). 

Plugging this into (IA.2j) with A so that Ai 7^ ±A 2 (which exists if N ^>d 1) 
gives 

(A.5) (aiAl ) 2 _(a 2 A 2 ) 2 = ^(/5 2 -a 2 ) 

and replacing A = (Ai, A 2 , . . . ) with (A 2 , Ai, . . . ) gives 

(A.6) ( ai A 2 ) 2 -(a 2 A 1 ) 2 = ^(/? 2 -a 2 ). 

Comparing (jA.5|) with (|A.6h gives 

(aiAi) 2 - (a 2 A 2 ) 2 = ( ai A 2 ) 2 - (a 2 Ai) 2 

that is 

(A 2 -A 2 )(a 2 + a 2 ) = 0. 

Since we chose A 2 7^ ±Ai this gives ai = a 2 = 0, contradicting ai 7^ 0. Thus 
we are done with step 1. 
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Step 2: We take C = ±D and a = ±/3 in (jA.ip and wish to show that for 
almost all z G T d , the only solutions are a = and C = 0. If either {a = (3 
and C = D) or (a = -/? and C = — -D), then (jSZQ) gives e 27ri ^' A > = -1 
which is a measure zero condition. 

Otherwise, assume a = (3 and C = —D (the other case is treated simi- 
larly). Here we have 

(A.7) a + i(C, A) = - e 2ni{x ' z) (a - i{C, A)) . 

If a = and C / then there is some A G A so that (C, A) / and (|A7fjl 
forces e 27ri ( 2 ' A ^ = 1, that is z lies on one of the hyperplanes 

Ua £ a{^ : (A, z) = mod 1}, 

which is a measure zero condition. 

If a 0, we replace C by — — C and a by 1 and drop the negative sign. 
Taking the real part of (|A,7j) . we have 

1 + cos 2vr(A, z) = (C, A) sin 2vr(A, z) . 

We may assume that the sine on the RHS doesn't vanish, since sin27r(A, z) = 
is a measure zero condition. Therefore, we may divide to get 

. . 1 + cos 2ir(\,z) 

(C, A) = ; = COt 7T(A, Z) . 

x ' 1 sin2vr(A,z) x ' 1 

Now square and average the result over an orbit O C A of the group Wd 
of all permutations and sign changes of the coordinates. The LHS gives 

1 1 xeo 

by (12. 3p , which is independent of the orbit chosen. The RHS gives 

T^TT COt 2 7r(A, Z) = 1 + — 7- 7-77 

^ x 1 O ^ sinvr (\,z)) 2 

1 1 xeo xeo y s ' " 

that is we find 

r^iri 2 - 1 = — V 1 

d l 1 |0|^(sinvr(A,z)) 2 • 

Since l/(sin 7r(A, z)) 2 is even, we get the same term for A and —A and so 
we may replace the average over O by the average over 0/ ± where we have 
taken only one of A, — A. Thus 

(A ' 8) 7 |C|2 " 1 = W±\ JC ± (sinvr(A,,)) 2 " 

Assuming that M > \Wd\ = 2 d d\, we can find a different orbit C C A 
and then comparing with (|A,8p gives 

1^1 1^1 



10/ ±1 ^ (sin^(A,z)) 2 107 ±1 ^ (sin7r(A,z)) 2 

that is we have eliminated the variable C. 

We claim that (|A.9p forces the point z to lie on a measure zero subset of 
T d . Indeed, the functions involved are meromorphic in C d /Z, d and hence if 
(|A.9p does not hold for all z, it can only hold on a complex submanifold of 
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codimension (at least) one and in particular its real points will have codi- 
mension at least one in T d . But near the origin z = 0, each of the functions 
l/(sin7r(A, z)) 2 has singularities on the hyperplane {z, A) = and these hy- 
perplanes are distinct for A's which are not collinear (here the condition 
d > 2 comes in), as is the case for those appearing in (LQ1) . Thus these 
functions are linearly independent and so (|A.9P is not valid for all z. □ 
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